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ABSTRACT 

We present a study of the multi-wavelength properties, from the mid-infrared to the hard X-rays, of a sample of 255 spectroscop- 
ically identified X-ray selected Type-2 AGN from the XMM-COSMOS survey. Most of them are obscured the X-ray absorbing 
column density is determined by either X-ray spectral analyses (for the 45% of the sample), or from hardness ratios. Spectral Energy 
Distributions (SEDs) are computed for all sources in the sample. The average SEDs in the optical band is dominated by the host-galaxy 
light, especially at low X-ray luminosities and redshifts. There is also a trend between X-ray and mid-infrared luminosity: the AGN 
contribution in the infrared is higher at higher X-ray luminosities. We calculate bolometric luminosities, bolometric corrections, stel- 
lar masses and star formation rates (SFRs) for these sources using a multi-component modeling to properly disentangle the emission 
associated to stellar light from that due to black hole accretion. For 90% of the sample we also have the morphological classifications 
obtained with an upgraded version of the Zurich Estimator of Structural Types (ZEST+). We find that on average Type-2 AGN have 
lower bolometric corrections than Type- 1 AGN. Moreover, we confirm that the morphologies of AGN host-galaxies indicate that there 
is a preference for these Type-2 AGN to be hosted in bulge-dominated galaxies with stellar masses greater than 10 10 solar masses. 

Key words, galaxies: active - galaxies: evolution - quasars: general - methods: statistical 

1 . Introduction ture formation over cosmic time. These models invoke different 

mechanisms to fuel the central SMBHs and to build the host- 

The formation and growth of supermassive black holes ga i axv bulges, such as major/minor mergers of gal axies (e.g. , 

(SMBHs) and their host-galaxies are related processes. This Corbin 2000; Kauffmann & H aehnelt 2000; Smingel et al. 2005; 

is supported by various observational signatures: the SMBH Hopkins et al 2006)> sm0Qth accretion of cold gas from fila . 

mass correlates with the_mass of the bulge ofthe host- me ntarv structures (e.g.. iKeTel^ l2009l Ir5ekel^d]l2 009). 

galaxy d Magorrian et al. | [19981 Marconi & Hunt 2003|), with Qr accretion of recycled gas from dying star s (e.g.. ICiotti^ttll 

flievelocity dispersionof the bulge ( |Ferrarese & Merrill] j^jJJ Several works also consider radiative fee dback which 

|200fJ Tremaineetal.||2002|), andwith the luminosity of the can repr oduce two important phases of galaxy evolution, namely 

bulge dKormendy & Richstone| [1995b. From theoretical mod- an obscure d-cold-phase, when the bulk of star formation and 

els, AGN seem to be able to switch off cooling flows in black hole accretion occurSj and the following quieS cent hot 

clusters (Hoeft & Briiggen 2004) and star formation in galax- phase in which a ccret ion remains hig hly sub-Eddington and 

les d Somervilleetal. ||2008D, with the result that the SMBH unobscured (e.g., ISazonov et l2005r. IlT^o & Ciottil 1201 1). 

mass is related to the host-galaxy bulge mass (or vice-versa). In some of these modeK ±e obscured/unobscured AGN di- 

Feedback between an accreting SMBH and the host-galaxy chotomy isrnorerelated to two different phases of galaxy evo- 

may play an important role in galaxy formation and evolu- lution (Horjkins et al 2008) ^ rather than to an orie ntation effect 

tion. Understand the role of feedback is a demanding prob- ,j g un j^ ed mode i sc heme) 
lem for both observers and theorists. Semi-analytical models 

and hydrodynamical simulations have been developed to attempt The obscured/unobscured time dependent AGN dichotomy 

to link the formation and evolution of SMBHs to the struc- could be related to the bimodality in the rest-frame color dis- 

tribution of host- 



elisabeta.lusso2@unibo.it iNandra et alJ 1200 



-g alaxies f Rovilos & Georgantop oulosl 120071 ; 
1 iBrusa et all l2009t ISilverman et all f2009). 
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namely the red-sequence (or "red-cloud") and blue-cloud galax- 
ies. Broad-line AGN (if the morphology of the host-galaxy is 
available) are likely to be associated to galaxies belonging to 
the blue-cloud, while obscured objects to red passive galaxies. 
The green valley should be populated by transition objects. The 
picture above is probably a too crude approximation. Moreover, 
one should note that red sequence galaxies may well be pas- 
sively evolving galaxies without si g nificant star formati on (e.g., 
iRovilos & Georgantopoulosl 120071: iNandra et ap 120071) . rather 
than dusty starforming objects (e.g jBrusa et alJ l2009). 

Disentangling the contribution of the nuclear AGN from the 
host-galaxy properties in the broad band SED is fundamental 
to constrain the physical evolution of AGN and to place them 
into the context of galaxy evolution. In the standard picture the 
AGN energy output is powered by accretion onto SMBHs. The 
disk accretion emission is visible in the optical-UV as the blue- 
bump feature. The X-ray emission is believed to be due to a 
hot-electrons corona that surrounds the accretion disk, while 
the infrared emission is likely due to the presence of a dusty 
torus around the disk at few parsec from the center, which re- 
processes the n uclear radiation. According to the unified m odel 
of AGN (e.g., lAntonuccil 1 1993. lUrrv & Padovanjl 1 19951) . hot 
dust is located in the inner edge of the torus. However, recent 
studies predict and observe exceptions to the unified model. 
From the theoretical point of view, an alternative solution to 
the torus is the disk-win d scenario (e.g.. lEmmering et al.ll 1992t 
lElitzur & Shlosmanl f2006). From the observational side, AGN 



without any detectable hot dust em ission (e.g..lJiangetalJ2010h 
and weak infrared emission fe.g.. lHao et alJl2010h are predicted 
and observed. The vast majority of studies performed so far con- 
cern unobscured (Type-1) AGN for which their SE D is well 
known from low-z ( (z) ~ 0.206, se e Elvis et aill 19941) to high- 
z ((z) ~ 1.525, see lRichards et al.ll2.Q06l) . An obvious compli- 
cation in the study of their host-galaxy properties is that the 
emission of the central AGN outshines the galaxy light in UV, 
optical and infrared bands; therefore it is extremely difficult to 
derive constraints on the colors, stellar populations, and mor- 
phologies of the host. On the other hand, for obscured (Type- 
2) AGN the host-galaxy light is the dominant component in 
the optical/near-infrared SED, while it is difficult to recover 
the AGN intrinsic nuclear emission. The lack of a proper char- 
acterization of the nuclear componentof the SED of obscured 
Type-2 AGN is a major limitation. As a consequence, the re- 
lations between stellar masses, SFR, morphologies and accre- 
tion luminosity remain poorly known. Since the relative con- 
tribution in the SED of the different components (AGN/host- 
galaxy) varies with wavelength, a proper decomposition can be 
obtained by an SED-fitting approach, complemented by a mor- 
phological analysis. This will provide a robust estimate of the 
nuclear emission (bolometric luminosities and bolometric cor- 
rections, absorption column density distributions, etc) and its re- 
lation with the host-galaxy properties (mass, star formation rates, 
morphological classification). The AGN structure is reflected in 
the shape of the SED, specifically the big-blue bump and the 
infrared-bump are related to the accretion disk and the surround- 
ing torus, respectively. Therefore, a densely sampled SED over 
a broad wavelength interval is mandatory to extract useful in- 
formation from SED fitting procedures, allowing to tightly con- 
strain physical parameters from multi-component modeling and, 
in particular, to properly disentangle the emission associated to 
stellar light from that due to accretion. 

The combination of sensitive X-ray and mid-IR observa- 
tories allows us to model the obscuring gas that in Type-2 
AGN hides the nuclear region from the near-IR to the UV. 



As supported by previous investigations, the reprocessed IR 
emission cou l d be a good proxy of the intrinsic disk emission. 
iGandhi et al.l (|2009) confirm the correlation between the X-ray 
luminosity at [2-10]keV and the IR emiss ion for a sample of 
Seyfert galaxies (see also iLutz et al J 120 04). Their data are the 
best estimate of the nuclear (non stellar) IR flux in AGN to date. 
A highly significant correlation between Lp-io]keV and the intrin- 
sic nuclear IR luminosity at 12.3^m is obse rved in the high qual - 
ity near-IR and X-ray data discussed by Gand hi et al.l (|2009). 
This reinforces the idea that "uncontaminated" mid-IR contin- 
uum is an accurate proxy for the intrinsic AGN emission. 

In this work we present the largest study of the multi- 
wavelength properties of an X-ray selected sample of ob- 
scured AGN using the XMM-Newton wide field survey in 
the COSMOS fi eld (XMM-COSM OS) Following a similar ap - 
proach to that of lPozzi et ai] _2007) and iVasudevan et al.l d2010). 
we use the infrared emission to evaluate the nuclear bolometric 
luminosity from a multi-component fit. The paper is aimed at 
a detailed characterization of a large sample of obscured AGN 
over a wide range of frequencies. The SEDs, morphology of the 
host-galaxies, stellar masses, colors, bolometric luminosities and 
bolometric corrections for the sample of obscured AGN are pre- 
sented. 

This paper is organized as follows. In Sect. [2] we report the 
selection criteria for the sample used in this work. Section [3] 
presents the multi-wavelength data-set, while in Sect. [4] the 
method to compute average SED is described. Section [5] con- 
cerns the multi-component modeling used to disentangle the nu- 
clear emission from the stellar light. In Section [6] the method 
used to compute intrinsic bolometric luminosites and bolomet- 
ric corrections for Type-2 AGN is described, while in Sect. Qwe 
have applied the same method to a sample of Type-1 AGN. The 
discussion of our findings is given in Sect. [8] while in Sect. [9] we 
summarize the most important results. 

We adopted a flat model of the universe with a Hubble con- 
stant Hp = 70 km s " 1 Mpc" 1 , Q M = 0.27, O a = 1 - Q 
(iKomatsu et al.ll2009b . 



a; 



2. The Data Set 

The XMM-COSMOS catalog comprises 1822 point-like X-ray 
sources detected by XMM-Newton over an area of ~ 2 deg 2 . The 
total exposure time was ~ 1.5 Ms with a fairly homogeneous 
depth of ~ 50 ks over a large fraction of the area (|Hasinger et al.l 
120071 ICappelluti et alj[2009b . Following iBrusa et all d2010l) , we 
excluded from our analysis 25 sources which turned out to 
be a blend of two Chandra sources. This leads to a total of 
1797 X-ray selected point-like sources. We restricted the anal- 
ysis to 1078 X-ray sources detected in the [2-10] keV band 
at a flux larger than 2.5 X 10~ 15 erg s _1 cm~ 2 (see Table 2 in 
ICappelluti et al.ll2009l) . The objects for which no secure optical 
counterpart could be assigned are often affected by severe blend- 
ing problems, so that we consider in this analysis the 97 1 sources 
(hereafter 971-XMM) for which a secure optical c ounterpart can 
be associated (see discussion in lBrusa et al.ll2010lPl 

From the 971-XMM catalog we have selected 255 sources, 
which do not show broad (FWHM< 2000 km s ) emission lines 
in their optical spectra^ (hereafter we will refer to them as the 



1 The multi-wavelength XMM-COSMOS catalog can be retrieved 
from: http://www.mpe.mpg.de/XMMCosmos/xmm53_release/, version 
1 st April 20 10~ 

2 The origin of spectroscopic redshifts for the 255 sources is as 
follows: 11 objects from the SDSS archive, 2 from MMT observa- 
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255 spectro-z 




Fig.l. Plot of the [2 - 10]keV flux versus the total i* CFHT 
magnitude for our sample of 255 Type-2 AGN. The red circles 
represent sources with a de-absorbed 2-10 keV luminosity lower 
than 10 42 erg s _1 . The dashed lines represent a constant X-ray 
to optical flux ratio of Log (X/ O) = + 1 . 



Type-2 AGN sample): 223 are classified not-broad-line AGN, 
while 32 are absorption-line galaxies. Not-broad-line AGN are 
objects with unresolved, high-ionization emission lines, ex- 
hibiting line ratios indicating AGN activity, and, when high- 
ionization lines are not detected, or the observed spectral range 
does not allow to construct line diagnostics, objects without 
broad line in the optical spectra. Absorption-line galaxies are 
sources consistent with a typical galaxy spectrum showing only 
absorption line. 

In Figure Q] we plot the 2-10 keV X-ray flux as a func- 
tion of i* CFHT magnitude. The dashed lines limit the region 
typically occupied by AGN along the X-ray to optical flux ra- 
tio Log (X/O) = +10. Nine sources have a de-absorbed 2- 
10 keV luminosity lower than 10 42 erg s , the conventional 
threshold below which the X-ray sources that can plausibly 
be explained by moderate-strength starbursts, hot gas in ellip- 
tical gal axies, or other sources bes ides accretion onto a nuclear 
SMBH (Hornschem eier et al.l200ll) . The three sources inside the 
dashed lines have X-ray luminosities close to 10 42 erg s _1 , while 
six AGN (6/255, 2%) lie in the part of the diagram usually oc- 
cupied by star-forming galaxies, and have X-ray luminosities 
< 10 42 erg s _1 . Their inclusion in the analysis does not affect the 
main results. 

The Type-2 AGN sample used in our analysis comprises 255 
X-ray selected AGN, all of them with spectroscopic redshifts, 
spanning a wide range of redshifts (0.045 < z < 3.524) and X- 
ray luminosities (41.06 < Log Lp-iojkev ^ 45.0). The redshift 
distribution of the total sample and the distribution of the de- 
absorbed hard X-ray luminosities are presented in Figure|2] The 



mean redshift is (z) = 0.76, while the mean Log Lp-tojkeV is 
43.34 with a dispersion of 0.64. 

2. 1 . Correction for absorption for the X-ray luminosity 

For a sub-sample of 1 1 1 AGN we have a n estimate of the col - 
umn density A^h from spectral analysis (see Mai nieri et al.l2007l) . 
while for 144 AGN absorption is estimated from hardness ra- 
tios (HR; see iBrusa et al.1 12010). For 25 sources for which a 
column density estimate is not available from HR we consider 
the Galactic column density. Therefore, we can compute the 
de-absorbed X-ray luminosity at 0.5-2 keV (soft band) and 2- 
10 keV (hard band) for all sources in our sample. In Figure [2] 
(right panel) we show the distribution of column densities that 
ranges from 3 x 10 20 cirT 2 to 1.5 x 10 24 cirT 2 . The mean A^h 
value is 8.5 x 10 21 cirT 2 with a dispersion of 0.72 dex. The in- 
tegrated intrinsic un-absorbed luminosity is computed assuming 
a power-law spectrum with slope, F = 2 and T = 1 .7 for the 
0.5-2 keV and 2-10 keV bands, respectively. The average shift 
induced by the correction for absorption in the Type-2 sample is 
<ALogL [2 _io]kev> = 0.04 ±0.01. 



3. Rest-frame monochromatic fluxes and 
Spectral Energy Distributions 



tions (jPresco tt et al. 2006), 70 from the IMACS observation campaign 
dTrump et alj|2q07h . 156 from the zCOSMOS bright 20k sample (see 
iLillv et alj|2007h . 7 from the zCOSMOS faint catalog and 9 from the 
Keck/DEDVIOS campaign. 
3 Log (X/O) = Log f x + f/2.5 + 5.6. 



We used the catalog bv lBrusa et al .1 (120101) which includes multi- 
wavelength data from mid-infrare d to hard X-ray s: MIPS 160 
yum, 70 //m an d 24 um GQ3 data £"e Floc'h et al.ll2009l) . IRAC 
flux densities (Sanders et alJ l2007h7 u* and i* CFHT bands and 
near-infrared K s -b and data dMcCracken et al.ll2008l) . J UKIRT 
dCapak et al 120081). HST /ACS F814W imaging of the COSMOS 
field dKoekemoer et al.l 120071). op tical multiband photometry 
(SDSS, Subaru, Capak et al. 2007) and ne ar- and far-ultraviolet 
bands with GALEX (IZamoiski et al.ll2007l) . 

The number of X-ray sources detected at 160 fxm and 70 fim 
is 18 and 42, respectively. For the undetected sources in these 
bands we consider 5<x upper limits of 65 mJy and 8.5 mJy for 
160 fim and 70//m, respectively. At 24 pmi the number of de- 
tected sources is 237. For the 18 undetected sources at 24fim, 
we consider 5<x upper limits of 80^/Jy. All 255 sources are de- 
tected in the infrared in all IRAC bands, and only a few objects 
were not detected in the optical and near-IR bands: we have only 
8 upper limits in the z + band; 1 upper limit in the Bj and i* 
bands; 2 upper limits in the u* band; 4 upper limits in the Ks 
CFHT band and 2 in the / UKIRT band. The observations in 
the various bands are not simultaneous, as they span a time in- 
terval of about 5 years: 2001 (SDSS), 2004 (Subaru and CFHT) 
and 2006 (IRAC). Variability for absorbed sources is likely to 
be a negligible effect, but, in order to further reduce it, we se- 
lected the bands closest in time to the IRAC observations (i.e., 
we excluded SDSS data, that in any case are less deep than other 
data available in similar bands). GALEX bands are not taken 
into account because, given the large aperture they can include 
light from close companions. All the data for the SED compu- 
tation were shifted to the rest frame, so that no K-corrections 
were needed. Galactic reddening has been taken into account: 
we used the selective a ttenuation of the st ellar continuum k(A) 
taken from Table 1 1 of Capak et al.l (120 07). Galactic extinction 
is estimated from ISchlegel et a" dl998l) ~for each object in the 
971-XMM catalog. Count rates in the 0.5-2 keV and 2-10 keV 
are converted into monochromatic X-ray fluxes in the observed 
frame at 1 and 4 keV, respective ly, using a Galactic col umn den- 
sity N H = 2.5 x 10 20 ct7j- 2 (see Capp elluti et"aDl2009t) . and as- 
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Fig. 2. Left panel: Redshift distribution of the 255 Type-2 AGN considered in this work. Center panel: Intrinsic hard X-ray lumi- 
nosity distribution of the 255 Type-2 AGN considered in this work. Right Panel: Column density distribution of the 255 Type-2 
AGN (black histogram), of the 111 Type-2 AGN with an Nh estimate from spectral analysis (grey filled histogram), and of the 144 
Type-2 AGN with an Nh estimate from hardness ratios (orange hatched histogram). 




Log(i/) [Hz] 

Fig. 3. Mean (orange crosses) and median (red points) SED from 
the total sample of 255 Type-2 AGN. The blue points represent 
the rest-frame data, from infrared to X-ray, used to construct 
the average SED, while the black lines represent the interpolated 
SED. 



suming a photon index F v = 2 and Y x = 1.7, for the soft and 
hard band, respectively. We do not correct these X-ray fluxes 
for the intrinsic column density. All sources are detected in the 
2-10 keV band by definition of the sample, while in the soft band 
we have 70 upper limits. 

4. Average SED 

We have computed the individual rest-frame SEDs for all 
sources in the sample, following the same approach as in L10, 
and we have normalized all of them at 1/im. After this nor- 
malization we divided the frequency interval from Log v^,, to 
Log v max using a fixed step ALog v. The minimum and max- 
imum frequency depends on both the data and the redshift of 
the source considered to compute the SED, in our case we 



have used Log v m j n = 12 Hz and Log v max = 20 Hz, with 
a ALog v = 0.02. We averaged data in each given interval 
Log v/ < Log v < Log v; + i . The mean and median SEDs are ob- 
tained by taking the arithmetic mean and the median of logarith- 
mic luminosities, Log L, in each bin. It is important to note that 
sources at different redshift contribute to the same bin. Because 
of the relatively wide range of redshifts, the lowest and the high- 
est frequency bins are populated by a variable number of points. 
This effect may introduce relatively high fluctuations in the av- 
erage luminosity in those bins. In order to minimize these ef- 
fects we select a minimum number of 200 SEDs in each bin (this 
number depends on the total number of sources in the sample). 
Then, we select the mean reference frequency, Log v, of the bin 
and use a binary- search algorithm to find all luminosities that 
correspond at Log v (if a source does not have a frequency that 
correspond to Log v, we choose the luminosity with the closer 
frequency to Log v). Finally, all adjacent luminosities in each 
bin are then connected to compute the final mean and median 
SED. In Figure [3] the resulting mean and median SEDs are re- 
ported with orange crosses and red points, respectively. The data 
points are reported in order to show the dispersion with respect 
to ljum. The average SED is characterized by a flat X-ray slope, 
(F = 1.12), while in the optical-UV the observed emission ap- 
pears to be consistent with the host-galaxy. The flattening of the 
X-ray slope is likely due to the fact that we do not correct for 
the intrinsic absorption the fluxes at 1 and 4 keV. The average 
SED in the mid-infrared is probably a combination of dust emis- 
sion from star-forming region and AGN emission reprocessed 
by the dust. Before trying to deconvolve each source using an 
SED-fitting code, we binned the total sample in X-ray and in- 
frared luminosites and redshift. We used the luminosity at 4 keV 
and 8jum to divide the total sample in 6 bins with the same num- 
ber of sources in each bin. The wavelength of the luminosity 
used to bin the sample is chosen to minimize the host-galaxy 
contribution. In the three panels of Figure [4] the resulting mean 
SEDs are shown. The redshift trend is directly connected to the 
luminosity trend, since at higher redshifts we are looking at the 
more luminous sources. The shapes of the average SEDs in the 
optical bands are approximately the same in all luminosity and 
redshift bins. As expected, there is a stronger host-galaxy con- 
tribution at lower luminosity/redshift bins, where the average 
SEDs have a typical galaxy shape. Moreover, there is a trend 
between X-ray and mid-infrared luminosity: the contribution in 
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Fig. 4. Average SEDs in the rest-frame Log (vL v ) - Log v plane. Right panel: Mean SEDs computed binning in X-ray luminosity 
at 4 keV. Center panel: Mean SEDs computed binning in infrared luminosity at 8/im. Left Panel: Mean SEDs computed binning in 
redshift. The color code refers to the different bins as labeled. 



the infrared is higher at higher X-ray luminosities. This effect 
is already known for both Type-1 and Type-2 AGN using the 
intrinsic (non-stellar) em ission from the AGN (e.g., iLutz et all 
12004 iGandhi et al1l2009h . The observed average SED for our 
sample is a combination of both host-galaxy and AGN light, but 
the change in the average shape in the mid-infrared as a function 
of the X-ray luminosity suggests that most luminous sources are 
probably dominated by the AGN emission at those wavelengths. 

5. SED-fitting 

The purpose of performing SED fitting is to properly disentangle 
the emission associated to stellar light from that due to accretion 
and constrain physical parameters. Since the relative contribu- 
tion of the different components varies with wavelength, a proper 
decomposition can be obtained through an SED-fitting approach 
providing a robust estimate of the nuclear emission (bolometric 
luminosities and bolometric corrections) and its relation with the 
host-galaxy properties (mass, star formation rates, morphologi- 
cal classification). A well sampled SED is mandatory; in par- 
ticular, far-infrared observations are fundamental to sample the 
star-formation activity, while mid-infrared observations are nec- 
essary to sample the region of the SED where most of the bolo- 
metric luminosity of obscured AGN is expected to be re-emitted. 
In Fig. the broad-band SEDs of four XMM-Newton Type-2 
AGN are plotted as examples. The lower two panels are repre- 
sentative of a full SED with all detections from the far-infrared 
to the optical. Unfortunately, there are very few detections at 160 
and 70yum, so that the more representative situation is shown in 
the upper left panel of Fig. [5] The three components adopted 
in the SED-fitting code, starburst, AGN torus and host-galaxy 
templates, are shown as a long-dashed line, solid line and dotted 
line, respectively. All the templates used in the SED-fitting code 
will be described in the following Sections. The red line repre- 
sents the best-fit, while the black points represent the photomet- 
ric data used in the code, from low to high frequency: MIPS- 
Spitzer (\6Qpm, 70pm and 24pm if available), 4 IRAC bands, K 
CFHT, J UKIRT, optical Subaru and CFHT bands. 

The observed data points from infrared to optical are fitted 
with a combination of various SED templates (see Sect. 15. U us- 
ing a standard^ 2 minimization procedure 

2 X" 1 F b s j — A X F gal J — B X F a gnJ — C X F [ r [ 



where F h s j and <x, are the monochromatic observed flux and 
its error in the band i; F ga ij, F ngn j and F irji are the monochro- 
matic template fluxes for the host-galaxy, the AGN and the star- 
burst component, respectively; A, B and C are the normaliza- 
tion constants for the host-galaxy, AGN and starbust compo- 
nent, respectively. The starburst component is used only when 
the source is detected at 160/im and 70pm. Otherwise, a two 
components SED-fit is used. Sixteen is the maximum number of 
bands adopted in the SED-fitting (only detection are considered), 
namely: 160fim, 70pm, 24pm, 8.0pm, 5.8pm, 4.5pm, 3.6pm, K$, 
J, z + , i", r + , g + , Vj, Bj and u*. 



5.1. Template libraries 
5.1.1. Optical template library 

We used a set of 75 galaxy templates built from the 
iBruzual & Cha riot (2003, BC03 hereafter) code for spectral syn- 
thesis models, using the version wi th the "Padov a 1994" tracks, 
solar metallicity and Chabrier IMF (Chabrier 2003). For the pur- 
poses of this analysis a set of galaxy templates representative of 
the entire galaxy population from ellipticals to starbursts is se- 
lected. To this aim, 10 exponentially decaying star formation his- 
tories with characteristic times ranging from t = 0.1 to 30Gyr 
and a model with constant star formation are included. For each 
SFH, a subsample of ages available in BC03 models is selected, 
to avoid both degeneracy among parameters and speed up the 
computation. In particular, early-type galaxies, characterised by 
a small amount of ongoing star formation, are represented by 
models with values of t smaller than 1 Gyr and ages larger than 
2Gyr, whereas more actively star forming galaxies are repre- 
sented by models with longer values of t and a wider range of 
ages from 0.1 to 10 Gyr. An additional constraint on the age is 
that, for each sources, the age has to be smaller than the age of 
the Universe at the red shift of the source. E ach template is red- 
dened according to the lCalzetti et a l. (2000) reddening law. The 
input value is E(B - V), corresponding to a dust-screen model, 
with F (A) = Fi(A)lQ-° AE{B - v>kw , where F and F t are the ob- 
served and the intrinsic fluxes, respectively. The extinction at a 
wavelength A is related to the colour excess E(B - V) and to the 
reddening curve k(A) by A A = k(A)E(B -V) = k(A)A v /R v , with 
Rv = 4.05 for the Calzetti's law. The E(B - V) values range 
between and 1 with a step of 0.05. 
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5.1 .2. AGN template library 

The nuclear SED templates are taken from ISilva et al.l (|2004). 
They were constructed from a large sample of Seyfert galaxies 
selected from the literature for which clear signatures of non- 
stellar nuclear emission were detected in the near-IR and mid-IR. 
After a proper subtraction of the stellar contribution, the nuclear 
infrared data were interpolated with a radiative transfer code for 
dust heated by a nuclear source with a typical AGN spectrum, 
and including different geometries, dust distribution, variation 
of the radii, density and dust grain sizes to account for possible 
deviati ons from a standard ISM extinction curve (se e for more 
details iGranato & Danesdll994t iMaiolino et al.ll2001l) . 

The infrared SEDs were then normalized by the intrinsic, 
unabsorbed X-ray flux in the 2-10 keV band, and are divided 
into 4 intervals of absorption: A^h < 10 22 cirT 2 for Syl, while 
10 22 < N H < 10 23 cirT 2 , 10 23 < N H < 10 24 cm' 2 and Na > 
10 24 cirT 2 for Sy2 (see Fig. 1 in ISilva et al.ll2004 . The main 
differences between the SEDs of Syls and Sy2s with 10 22 < 
A^h < 10 23 cm~ 2 are the absorption in the near-IR at about A < 
2yum and the silicate absorption at A = 9.7//m, which are present 
in the Sy2 template. The shape of the SED in the mid-infrared 
with 10 23 < jV h < 10 24 cirT 2 is quite similar to that with 10 22 < 
N H < 10 23 cm" 2 . The Compton-thick SED (N H > 10 24 cirT 2 ) 
shows conspicuous absorption also at A ~ 1 .3/itn. If a source has 
the A^h value available, this is used as a prior in the selection of 
the best-fit AGN template. 



5.1.3. Starburst template library 

We us ed two different starbur st te mplate libraries for th e SED- 
fitting: ICharv & Elbazl (l200ll) and lDale & Heloul (I2002I) . These 
template libraries represent a wide range of SED shapes and lu- 
minosities and are widely used in the literature. Here, we briefly 
describe how each of these libraries was derived and discuss the 
main differences between them. 

The ICharv & Elbazl (l200ll) template library consists of 105 
templates based on the SEDs of four prototypical starburst galax- 
ies (Arp220 (ULIRG); NGC 6090 (LIRG); M82 (starburst); and 
M5 1 (normal star -forming galaxy)). They were derived using the 
ISilva et ail d 1998b models with the mid-infrared region replaced 
with ISOCAM observations between 3 and 18/mi (verifying that 
the observed values of these four galaxies were reproduced by 
the templates). These templates were then divided into two por- 
tions (4-20/im and 20-1000/im) and interpolated between the 
four to gener ate a set of librar ies of varying shapes and lumi- 
nosities. The iDale et all d2001l) templates are also included in 
this set to extend the range of shapes. 

The lDaie& Heloul d2002l) templates are updated versions of 
the lDale et al.l(l200ll) templates. This model involves three com- 
ponents, large dust grains in thermal equilibrium, small grains 
semistochastically heated, and stochastically heated PAHs. They 
are based on IRAS/ISO observations of 69 norm al star-forming 
galaxies in the wavelength range 3-100yum. IDale & Heloul 
(2002) improved upon these models at longer wavelengths using 
SCUBA observatio ns of 114 galaxies from the Bright Galaxy 
Sample (BGS, see ISoifer et al.l Il989l) . 228 galaxies observed 
with ISOLWS (52-170/im; Brauher 2002), and 170//m obser- 
vations for 115 gala xies from the ISOPHOT Serendipity Survey 
(IStickel et al.1 12000). All together, these 64 templates span the 
IR luminosity range 10 8 - 10 I2 L Q . The total infrared template 
sample used in our analysis is composed of 168 templates. 



6. Bolometric luminosities and bolometric 
corrections 

The nuclear bolometric luminosities and bolometric corrections 
are estim ated, using an approach similar to Pozz i et al. (2007, 
see also Vasu devan et al.ll2010fc iPozzi et al.ll2010l) . whereas the 
infrared luminosity is used as a proxy of the intrin sic nuclear 
lumino sity. The appropriate nuclear template from ISilva et aT] 
(l2004h is selected based on the absorbing column density A^h, 
when available, or from the best-fit nuclear infrared template. In 
order to compute the hard X-ray bolometric correction we used 
the standard definition 



M2-10]keV 

where the Lp-iojkev is the intrinsic X-ray luminosity and the 
bolometric luminosity is computed as the sum of the total in- 
frared and X-ray luminosity 

Lboi = L\r + L x . (3) 

After performing the SED-fitting, only the nuclear component of 
the best-fit is integrated. Hence, the total IR luminosity Lir is ob- 
tained integrating the nuclear template between 1 and 1000/im. 
To convert this IR luminosity into the nuclear accretion disk 
luminosity, we applied the correction fact ors to account fo r 
the torus geometry and the anisotropy (see IPozzi et alj|2007l) . 
The first correction is parameterized by the covering factor /. 
The covering factor is related to the geometry of the torus that 
obscures the accretion disk emission in the optical-UV along 
the line of sight, and its value is estimated from the ratio of 
obscured/unobscure d quasars foun d by the X-ray background 
synthesis models dGilli et al.l 120071) . This correction factor is 
~ 1.5. This value correspond to a typical covering factor of 
/ ~ 0.67, consistent with the results based on clumpy torus mod- 
els dNenkova et al.ll2008l) . 

The anisotropy factor is defined as the ratio of the luminosity 
of face-on versus edge-on AGN, where the obscuratio n is a func- 
tion o f the column density A^h- Therefore, SEDs in Sil va et al.l 
(2004) have been integrated in the l-30//m range, after normal- 
izing these SEDs to the same luminosity in the 30-100/im range. 
The derived anisotropy values are 1.2-1.3 for 1 22 < A^h < 10 24 
and 3- 4 for A^h > 10 24 . The same values as in IVasudevan et al] 
(l2010h are adopted: 1.3 for 10 22 < N H < 10 24 and 3.5 for 
N H > 10 24 . 

The total X-ray luminosity Lx is estimated integrating in the 
0.5-100 keV range the X-ray SED. We have interpolated the 
de-absorbed soft and hard X-ray luminosities. Since we are in- 
tegrating at the rest-frame frequencies, the X-ray SED is ex- 
trapolated at higher and lower energies using the estimated X- 
ray slope, and in troducing an exponential cut-off at 200 keV 
dGilli et al.ll2007l see also Sect. 3.1 in L10). 

7. Robustness of the method 

The robustness of the method used to estimate nuclear bolomet- 
ric luminosities and bolometric corrections from SED-fitting, for 
the sample of Type-2 AGN, has been tested against the updated 
soft X-ray selected sample of Type-1 AGN discussed in L10. 
The Type-1 AGN sample in the L10 work was composed of 361 
spectroscopically classified broad-line AGN. The recent work by 
Brusa et al. (2010) has updated the spectroscopic classification 
and increased the number of Type-1 AGN with spectroscopic 
redshift, so that the final sample is composed of 395 Type-1 
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Fig. 5. Examples of SED decompositions. Black circles are the observed photometry in the rest-frame (from the far-infrared to the 
optical-UV). The long-dashed, solid and dotted lines correspond respectively to the starburst, AGN and host-galaxy templates found 
as the best fit solution. The red line represents the best-fit SED. The stellar mass and the SFR derived from the galaxy template are 
reported. 



AGN in the redshift range 0.103 < z < 4.255 with X-ray lu- 
minosities 42.20 < Log Lp-iojkev ^ 45.23. We have computed 
bolometric and X-ray luminosities, and bolometric corrections 
using the same approach as in L10 for the Type-1 sample: bolo- 
metric luminosites are computed by integrating the rest-frame 
SEDs from lpm up to the UV-bump. In order to compare these 
estimates with the results from the SED-fitting code, we have 
applied to the same sample the method described in Sect. [5] and 
[6] to estimate bolometric parameters. To be consistent with the 
selection criteria of the sample discussed in this paper, we have 
considered only AGN with X-ray detection in the hard band, re- 
moving from the main sample 87 Type-1 AGN with an upper 
limit at 2-10 keV. Moreover, for 2 sources the best-fit does not 
consider an AGN component, so we cannot compute the bolo- 
metric luminosities for them. The final test sample is composed 
of 306 Type-1 AGN in the redshift range 0.103 < z < 3.626 and 
X-ray luminosities 42.20 < Log Lp-iojkeV 2 45.04. 

In order to s elect the appropriate nuclear template from 
ISilva et al.l d2004l) . we consider the SED for Syl AGN (no cor- 
rection for anisotropy is necessary in this case) and the Sy2 SED 
with 10 22 < N H < 10 23 for 20 AGN that have N H in this range. 

We present the comparison between the values of L\, Q \ and 
£boi from L10 and this work in Fig. [6] The outlier in the bottom 
side of the plot, XID=357 at redshift 2.151 has Log k bo] = 1.95 



from L10 and Log k bo \ = 1-04 using the new approach, and 
presents large error bars in the 24fj.m detection, so that the total 
bolometric luminosity, computed using the infrared luminosity, 
is probably underestimated. The outlier in the right end of the 
distribution, XID=51 14 at redshift 0.212 (Log k bol = 2.82 from 
L10 and Log k bo i - 1.95 using the new approach) has detections 
at 160, 70 and 24/vm, Log A^h = 22.68 and Log Lp-iojkev = 
42.89. Probably this source is a star-forming galaxy, so that us- 
ing the L10 approach we included stellar emission in the esti- 
mate of the nuclear bolometric luminosity, thus overestimating 
the nuclear bolometric luminosity and, therefore, the bolometric 
correction. The last notable outlier in the top/left side of the dis- 
tribution, XID=2152 atredshift 0.627 (Log £ bol = 1.35fromL10 
and Log fcboi = 2.27 using the new approach) presents a signif- 
icant host-galaxy contribution in the optical-UV and, therefore, 
the bolometric luminosity is likely to be underestimated in the 
L10 approach. 

Although the two methods are very different, the bolometric 
luminosity estimates agree remarkably well, with a lcr disper- 
sion of 0.20 dex after performing a 3.5<x clipping method in or- 
der to avoid outliers. Bolometric luminosities from SED-fitting 
are on average slightly larger than those computed integrating 
the rest-frame SED from \\xm to the X-ray (see the lower left 
side in Fig. [6]). This effect is also present in the Vasudevan et al. 
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Fig. 6. Upper panel: Comparison between the values of bolometric luminosity and bolometric correction from data presented in 
L10 and from this work. The three red triangles mark the outliers discussed in Sect. [7] Lower panel: Distribution of the differences 
between the values of bolometric luminosity and bolometric correction from data presented in L10 and from this work. 



(2010) work. A possible explanation is that SED-fitting underes- 
timates the host-galaxy contribution, or that the anisotropy and 
geometry corrections are too large for some objects. The agree- 
ment between the two methods is overall quite satisfactory and 
in the following we will discuss our findings for the Type-2 sam- 
ple. 

8. Results and discussion 

8.1. Bolometric correction and luminosites for Type-2 AGN 

Bolometric luminosities and bolometric corrections have been 
computed for the Type-2 AGN sample. Intrinsic soft and hard 
X-ray luminosities are estimated as described in Sect. I2.ll For 
15 sources we do not have an estimate of the AGN compo- 
nent from the SED-fitting, and we cannot compute the bolomet- 
ric luminosity for them. In Fig. [7] the distribution of the bolo- 
metric correction for the 240 Type-2 AGN sample and for the 
306 Type-1 AGN (febol f° r Type-1 AGN are computed using 
the SED-fitting code) are presented. Figure [8] shows bolomet- 
ric corrections for both the Type-1 and the Type-2 AGN sam- 
ples as a function of the hard X-ray luminosity. For both sam- 
ples, bolometric parameters are estimated from the SED-fitting 
as discussed in Sect. [5] The green and orange curves represent 
th e bolometric correctio ns a nd their 1 cr dispers ion as derived 
bv lHopkins et alj(l2007l) and lMarconi et al.l (120041) . respectively. 
Type-2 AGN have, on average, smaller bolometric corrections 
than Type-1 AGN at comparable hard X-ray luminosity. For 
example, at 43.30 < Log Lp-tojkev ^ 44.30 (vertical lines in 
Fig. [D, where both AGN types are well represented, the me- 




0.5 1 1.5 2 2.5 



Lo § k bol 

Fig. 7. Distribution of the bolometric correction for the 240 
Type-2 AGN sample (red hatched histogram) and for the 306 
Type-1 AGN (blue hatched histogram). 



dian bolometric correction for the Type-2 AGN (134 objects) is 
(^boi) ~ 13 + 1, to be compared with a median bolometric correc- 
tion (fc bol ) -23+1 for the Type-1 AGN (167 objects). The two 
averages are statistically different at the ~ 7 cr level and this is 
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Fig. 8. Hard X-ray bolometric correction against the intrinsic 
2-10 keV luminosity for 240 Type-2 AGN with AGN best-fit 
(red data). The crosses represent the bolometric correction for 
306 Type-1 AGN, computed with the approach described in 
Sect. [6] The green and blue lines represe nt the bolometric cor- 
recti on and the lcr dispe rsion obtained by Hopk ins et al.l (120071) 
and lMarconi et alJ (l2004h . respectively. The red points and open 
squares represent the 111 Type-2 AGN with Nh from spectral 
analyses and the 144 Type-2 AGN with Nh from HR, respec- 
tively. 



consistent with the results in lVasudevan et all (1201 Oh . The mean 
£p-i0]kev for the Type-1 and Type-2 AGN within this luminosity 
range differs by a factor 1.8, and this could in principle explain at 
least part of the difference in the average bolometric corrections 
for the two samples of AGN. However, the significance of the 
difference is still present if we split this luminosity range in two 
equal Log Lp-io]keV bins and perform a Kolmogorov-Smirnov 
test for the Type-1 and Type-2 AGN luminosity distributions in 

eac h bin. 

IVasudevan & Fabianl (l2009h and iLusso et al ] (1201 Oh have 
shown that hard X-ray bolometric corrections are correlated 
with the Eddington ratios (Amd = ^boi/^Edd) fo r Type-1 AGN 
(see also lMarconi et al.ll2004trKellv et al.ll2008l) . The /fc bol - A Edd 
relation suggests that there is a connection between the broad- 
band emission, mostly in the optical-UV, and the Eddington 
ratio, which is directly linked to the ratio between mass ac- 
cretion rate and Eddington accretion rate. A high AEdd corre- 
sponds to an enhanced optical-UV emission, which means a 
prominent big-blue bump and therefore a higher k\, \. The dif- 
ference between the average bolometric corrections for Type- 
1 and Type-2 AGN could be due to lower mass accretion rates 
in Type-2 AGN, assuming the same b lack hole mass dis tribu- 
tion for the two AGN populations (see lTrump et al.ll201 lh . The 
current theoretical framework of AGN/host-galaxy co-evolution 
predicts that obscured AGN are highly accreting objects and 
their black hole is rapidly growing. Howev er, we note that this 
is true for z — 1 - 3 (see Hasinger 2008), while the majority 
of Type-2 AGN in the present sample are relatively luminous 
(42.37 < Log Lboi ^ 46.80), and at moderately low redshift 
(0 < z < 1). Therefore, the Type-2 AGN sample is likely to 
represent a later stage in AGN evolution history. 



8.2. Infrared emission: indication of AGN activity 

The re-processed infrared emission can be used as a proxy of the 
average disc emission, since the timescale for transfer of energy 
from the disk to the outer edge of the torus into infrared emis- 
sion is of the order of several years in standard AGN picture; 
whereas optical, UV and X-ray variability in AGN is known to 
occur on shorter timescales. The correlation between the 2-10 
keV X-ray emission and IR emission at 12.3yum for a sample of 
Seyfert nuclei has been discussed in Gandhi et al. (2009), and it 
could be used to estimate the intrinsic AGN power. Using X-ray 
data from the literature and new IR data from the Very Large 
Telescope's Imager and Spectrometer for mid-Infrared (VISIR), 
taken specifically for addressing the issue of nuclear emission 
in local Seyferts, they found a tight correlation between intrin- 
sic, uncontaminated IR luminosity and X-ray luminosity in the 
2-10 keV range 



Log 



^12.3 fim 

~iaW 



(0. 1 9 + 0.05) + ( 1 . 1 1 + 0.07)Log 



^2-10]keV 

10 43 ' 



(4) 



The relation is characterized by a small scatter with a standard 
deviation of 0.23 dex. The expected nuclear mid-infrared lumi- 
nosity is computed from Eq. © using the estimate of the intrin- 
sic unabsorbed X-ray luminosity. From the observed rest-frame 
SED (AGN+host-galaxy) the luminosity at 12.3 //m is com- 
puted. A comparison of the total observed luminosity at 12.3 //m 
and that predicted by Eq. I© is plotted in Figure [9] for four rep- 
resentative sources. In Figure [10] the distribution of the ratio 

r = Log (L12.3 ^ m ,obs/ii2.3 ^.predicted) is plotted for the Type-2 
AGN sample. The distribution of the ratio r has a mean which is 
shifted from zero by ~ 0.2. However, if we consider a gaussian 
distribution cente red at r = w i th cr — 0.23, i.e. the same disper- 
sion observed bv lGandhi et al.l (|2009) in their local sample, the 
majority of the objects are found within 2<x of the r distribution. 
The tail outside 2<x and extending to high r includes 73 sources 
(with r > 0.5) for which the predicted mid-infrared luminosity is 
significantly lower than observed. The hard X-ray luminosities 
of these 73 AGN are mainly in range Log Lp-iojkev ~ 42 - 44, 
where the local correlation is well-sampled. There are two pos- 
sible explanations for a significant (~ 30% of the objects) tail 
toward high-r values: either the Gandhi relation, which was de- 
rived for a sample of local Seyfert galaxies, cannot be extended 
to all the sources in our sample or the SED-fitting procedure may 
overestimate, in a fraction of these objects, the nuclear contribu- 
tion. In order to study the properties of these outliers, bolometric 
corrections, morphologies, stellar masses and SFR are discussed 
in following. We call "low-r" AGN all sources within 2cr of the r 
distribution, while the "high-r" AGN sample is populated by the 
sources deviating more than 2cr (see Fig.[9]for some examples). 

A clear separation in bolometric corrections for these two 
sub-samples is found. This is shown in Figure[TT]in which bolo- 
metric corrections are plotted as a function of the 2-10 keV lumi- 
nosity. At a given hard X-ray luminosity (43 < Log Lp-iojkev ^ 
44) the low-r sample has a median bolometric correction of 
(^boi) ~ 11 ± 1 (110 objects), to be compared with a median 
bolometric correction for the high-r sample of (k\, \) ~ 26 ± 3 
(44 objects). The two median values for k\, Q i are statistically dif- 
ferent at the ~ 5cr level. 

Furthermore, in the high-r sample 24 Type-2 AGN out of 73 
have a detection at 70/mi (~ 33%, significantly higher than those 
for the low-r sample, ~ 4%) and 9 of these 24 have also a de- 
tection at 160/mi (~ 12% considering the total high-r sample, 
and only ~ 1% for the low-r sample). This denotes that the dif- 
ference in the average bolometric corrections between the low-r 
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Fig. 9. Examples of SED decompositions. Black circles are the observed photometry in the rest-frame (from the far-infrared to the 
optical-UV). The long-dashed, solid and dotted lines correspond respectively to the starburst, AGN and host-galaxy templates found 
as the best fit solution. The red line represents the best-fit SED. The stellar mass and the SFR derived from the galaxy template are 
reported. The green point represents the nuclear mid-infrared luminosity using Eq. ©, while the cross represents the total observed 
luminosity at 12.3 /mi computed from the rest-frame SED. XID=19 and 81 are examples of low-r AGN, while XID=172 and 117 
represent high-r AGN. 



and high-r samples is probably due to the fact that a significant 
fraction of the infrared emission is attributable to an incorrect 
modeling of the star-formation process, or the AGN contribution 
is somehow overestimated by the SED-fitting procedure. 

There is no significant difference in the average nuclear ab- 
sorption between the low-r and the high-r sample, while there 
is a possibly significant difference in SFR and stellar masses. 
The median stellar mass in the high-r sample is (Log M.) ~ 
1O.93M with a dispersion of 0.30, while for the low-r sam- 
ple is (Log M»> ~ 1O.76M with cr = 0.30. The two av- 
erages are statistically different at the ~ 3 cr level. The me- 
dian SFR, as derived from the SED-fit, for the high-r sample 
is (SFR) ~ 17 ± 3M /yrs with a cr = 0.30, while for the low-r 
sample is (SFR) ~ 3 + lM /yrs with a cr = 0.30 and the two 
averages are statistically different at the 4.4 cr level. 

Overall, the SED-fitting for the 73 Type-2 AGN is likely to 
overestimate the AGN emission in the infrared, which is proba- 
bly due to the infrared emission from star-forming regions. The 
average bolometric correction for Type-2 AGN, excluding these 
sources, would be even lower than what we have computed in 
the previous Section. This reinforces the idea of lower bolomet- 
ric corrections for Type-2 AGN with respect to Type-1 AGN. 



The low bolometric corrections for Type-2 AGN could be also 
explained if a fraction of the accretion disk bolometric output is 
not re-emitted in the mid-infrared, but rather dissipated (e.g., by 
AGN-feedback). This would not be accounted in the bolomet- 
ric luminosity, and could provide a plausible explanation for low 
&boi especially if the low-r sample is considered. At this stage, 
this is just a speculation, and more work is needed to verify this 
possibility. 

8.3. Host-galaxy properties: M„, SFR, colors and 
morphologies 

Galaxies show a colour bi-modality both i n the local Universe 
and at higher r edshift (up to z ~ 2; e.g., Stra teva et alJ [2001: 
iBell et al.l2004l) . This bi-modality (red-sequence and blue-cloud 
galaxies) has been interpreted as an evidence for a dicothomy 
in the ir star format ion and merging histor ies (e.g.. lMenci et al.l 
2005, but see also Cardamon e et al. l l2OT0l for an alternative ex- 
planation). Color-magnitude and color-mass diagrams (e.g., rest- 
frame (U - V) versus stellar mass) have been used as tools in 
galaxy evolution studies, and since many models invoke AGN 
feedback as an important player in such evolution, it is inter- 
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XID 


Redshift 


Log L[2-10]keV 


Log L bo] 




Log M, 


SFR 


Mu 


Mm 


M s 


Morphological class 3 






[ergs -1 ] 


[ergs" 1 ] 




[M G ] 


[M Q /yrs] 










67 


0.367 


42.79 


43.80 


10.24 


9.68 


1.86 


-18.63 


-19.80 


-20.72 





65 


0.979 


43.83 


44.85 


10.62 


10.19 


38.60 


-20.50 


-21.63 


-22.71 


7 


64 


0.686 


43.54 


44.52 


9.57 


10.45 


32.82 


-20.17 


-21.56 


-23.05 


10 


63 


0.355 


42.98 


44.28 


19.88 


10.79 


4.94 


-20.61 


-22.13 


-23.26 


2 


54 


0.350 


42.58 


43.38 


6.32 


11.18 


0.11 


-20.57 


-22.62 


-23.78 


12 


45 


0.121 


41.90 


42.91 


10.31 


9.39 


0.00 


-16.76 


-18.70 


-19.77 


1 


43 


1.162 


44.22 


45.28 


11.49 


11.30 


1.47 


-21.55 


-23.48 


-24.61 


2 


19 


0.659 


43.67 


44.76 


12.44 


11.07 


35.88 


-21.57 


-22.88 


-24.02 


1 


117 


0.936 


43.47 


45.59 


132.11 


11.24 


199.40 


-21.81 


-23.34 


-24.80 


23 


116 


0.874 


43.49 


44.50 


10.17 


10.56 


3.77 


-20.19 


-21.81 


-22.90 





112 


0.762 


43.65 


44.65 


9.97 


10.93 


0.62 


-20.00 


-22.13 


-23.52 


3 


104 


0.623 


44.08 


45.11 


10.75 


10.79 


0.45 


-20.25 


-22.18 


-23.30 


1 


101 


0.927 


43.70 


44.68 


9.56 


10.92 


0.61 


-21.19 


-22.93 


-23.78 





100 


0.270 


42.61 


43.54 


8.59 


10.21 


0.00 


-19.14 


-20.98 


-21.92 





99 


0.730 


43.54 


44.78 


17.47 


9.97 


116.26 


-21.61 


-22.26 


-23.15 


11 


85 


1.001 


43.46 


44.84 


23.99 


10.13 


38.52 


-20.15 


-21.34 


-22.56 


8 


81 


0.915 


44.11 


45.05 


8.65 


11.18 


0.00 


-20.35 


-22.59 


-24.07 


2 


70 


0.688 


44.00 


45.60 


39.74 


10.65 


542.60 


-20.87 


-22.32 


-24.31 


3 


152 


0.895 


43.75 


44.82 


11.72 


9.87 


92.88 


-21.06 


-21.82 


-22.85 


7 


150 


0.740 


43.29 


44.25 


9.27 


10.64 


0.55 


-19.38 


-21.28 


-22.50 


3 



Notes — This table is presented entirely in the electronic edition; a portion is shown here for guidance. 

a The morphological classification of the Type-2 AGN hosts is coded from to 23: = elliptical, 1 = SO; 2 = bulge-dominated; 3 = 
intermediate-bulge; 4 = disk-dominated; 5 = irregular; 6 = compact/irregular; 7 = compact; 8 = unresolved/compact; 9 = blended; 10 = 
bulge-dominated/close-companion; 11 = intermediate-bulge/close-companion; 12 = SO/close-companion; 23 = possible mergers. 




Fig. 10. Histogram of the ratio between the total observed lumi- 
nosity at 12.3 yum and the mid-infrared luminosity predicted by 
Eq. ©. The red curve represents a gaussian with mean equal to 
zero and standard deviation 0.23. The l<x and 2<x standard devi- 
ations of the correlation are also reported. 

esting to locate the hosts of Type-2 AGN in those diagrams. 
Using the galaxy component obtained from the best fit of the 
Type-2 AGN, it is possible to derive rest-frame colors for the 
host that, linked to the stellar mass and the morphology, can 
provide hints on AGN feedback. Several studies found that 
the hosts of obscured AGN tend to be redder than the over- 



Fig. 11. Hard X-ray bolometric correction against 2-10 keV lu- 
minosity for 240 Type-2 AGN with AGN best-fit. The 240 Type- 
2 sample is divided into subsamples: low-r AGN sample (black 
data) and high-r AGN sample (yellow data). The green and or- 
ange lines r epresent the bo l ometri c co rrection and lcr dispe rsion 
obtained by Hopki ns et ail d2007l) and lMarconi et all d2004l) . re- 
spectively. In the de-absorbed hard X-ray luminosity range high- 
lighted by the solid lines, we have 167 low-r and 73 high-r (in 
the infrared) sources. 



all galaxy populat ion in the rest-frame (U - V) color (e.g., 
Nand ra et al] 120071) . There are at least two possible and sig- 
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Fig. 12. The morphology distribution (using the ZEST+ code) 
of the 233 AGN host-galaxies on the (U - V) colour-mass di- 
agram. We also plotted the 22 sources without morphological 
information. The (U — V) color and stellar masses are computed 
using the SED-fitting code. We overplot the contours of about 
8700 galaxies in zCOSMOS (colours and stellar masses from 
Hyperz). The morphology classification is labeled as follow: el- 
liptical (Ell), SO, bulge-dominated galaxy (BD), intermediate- 
bulge galaxy (IB), disk-dominated galaxy (DD), irregular (Irr), 
Compact, possible mergers (PM) and unresolved compact (UC). 
The red dashed lin e represents the red sequence cut defined by 
iBorch et al.l ((2006), while the green short dashed line defines an 
approximate green valley region, both lines are calculated at red- 
shift ~ 0.76, which is the average redshift of the main Type-2 
sample. 



nificantly different interpretations for this observational result: 
the observed red colors are mainly due to dust extinction, so 
that a significant fraction of obscured AGN would live in mas- 
sive, dusty star-fo rming galaxies with red optical colors (e.g., 
iBrusa et aLll2009t); or red sources are linked with passive sys- 
tems (e.g.. iRov ilos & Georga ntopoulo s 2007; S chawinski et al.l 
2009; Cardam one et alj|2010l) . Therefore, accurate stellar mass 
and SFR estimates, together with detailed galaxy morphologies, 
are of particular importance to discriminate between the two al- 
ternative possibilities. 

The very high resolution and sensitivity of ACS-HST imag- 
ing in the COSMOS survey provides resolved morphologies 
for several hundreds of thousands galaxies with / acs < 24 (see 
Scarlata et al 2007 for details). Galaxy morphologies were ob- 
tained with an upgrade d version of the Z urich Estimator of 
Structural Types (ZEST; IScarlata et al.ll2007h . known as ZEST+ 
(Carollo et al. 2011, in prep). Relative to its predecessor, 
ZEST+ includes additional measurements of non-parametric 
morphological indices for characterising both structures and 
substructures. For consistency with the earlier versions, ZEST+ 
uses a Principal Component Analysis (PCA) scheme in the 
6-dimensional space of concentration, asymmetry, dumpiness, 
M20 (second-order moment of the brightest 20% of galaxy pix- 
els), Gini coefficient, and ellipticity. ZEST+ classifies galax- 
ies in seven morphological types located in specific regions 
of the 6-dimensional space: elliptical, SO, bulge-dominated 



disk, intermediate-bulge disk, disk-dominated, irregular, com- 
pact. The different types were then visually inspected. For 19 
objects ZEST+ is unable to give any information on morphol- 
ogy because these sources lie off the edge of the ACS tiles and 
4 sources are blended. As a result of the ZEST+ procedure and 
visual inspection of the other 233 galaxies in our sample, we find 
that 16 are ellipticals (Ell), 53 are SOs, 74 are bulge-dominated 
(BD) disks, 27 are intermediate-bulge (IB) disks, just 1 is disk- 
dominated (DD), 19 are irregular galaxies (Irr), 15 are compact 
galaxies (i.e. the structural parameters computed for these galax- 
ies from the HST-ACS images are highly affected by the in- 
strumental PSF) and 18 are unresolved compact galaxies (UC, 
i.e. essentially point-like sources). Ten galaxies show distortions 
and potential signatures of ongoing or recent mergers (PM). At 
the typical magnitudes of the objects in our sample, the ZEST+ 
classification is highly reliable for galaxies with redshift < 1. 
At higher redshifts, morphological k-correction and, to a lesser 
extent, resolution effects can adversely affect measurements of 
ZEST+ parameters (note that only 4% of the main Type-2 AGN 
sample have z > 1.5). However, broad morphological bins (e.g., 
early/late type galaxies) should be relatively robust. For high-z 
galaxies, resolution might also have an impact on the classifica- 
tion for mergers, and ACS images coul d not be deep enoug h 
to distinguish merger features (see also iMainieri et al. I l20!lh . 
Moreover, inc lination might also aff ect the morphological classi- 
fication (e.g jNair & Abrahaml2010l) . However, a detailed study 
of systematics and biases in the morphological classification is 
beyond the purposes of the present paper. In Table Q] we list the 
main properties of the sample. 

The rest-frame (U — V) color encompasses the 4000 A break, 
and it is particularly sentitive to age a nd metallicity variations 
of ste l lar populat i on in galaxies (e.g ., ISandage & Vi svanathan 
1 19781: iBell etalJ l2004t IBorch et al] 120061: ICardamone et all 
2010). In Fig.fTSIthe distribution of the rest-frame (U—V) colors, 
which are computed directly from the best-fit galaxy template, 
and stellar masses (from the SED-fitting code) are reported for 
the entire Type-2 AGN sample. In the same figure, the back- 
ground contours for a sample of ~ 8700 galaxies in zCOSMOS 
O'acs < 22.5, 240 Type-2 are detected in the z' acs band, 183/240 
Type-2 AGN (76%) have i acs < 22.5) are also plotted, where 
colours and stellar mas ses are computed using the Hyperz code 
dBolzonella et alJl2000h . 

AGN are known to reside in m assive galaxies (e.g., 
ISilverman" et all 120091: IBrusa et al.l F2009) and this is fully con- 
firmed by the present analysis. The morphologies of the host- 
galaxies and the stellar masses indicate that there is a prefer- 
ence for these Type-2 AGN to be hosted in bulge-dominated and 
SO galaxies (~ 50%) with stellar masses greater than 10 10 M o . 
This resul t is consistent with the prev ious studies on Type-2 
AGN bvlSilyerman et alj d2008l see also iKauffmann et al.l l2003b 
iBundv et al.ll2008l:ISchawinski et al.ll201 lb . 

It should be noted that no correction for the internal extinc- 
tion has been applied to the (U — V) colors of both background 
galaxies in zCOSMOS and Ty pe-2 AGN hosts. This correction 
could be important as sh own in lCowie & Ba rger (2008) (see also 
Cardamone et a 1. 2010). In that work star formation and galac- 
tic stellar mass assembly are analyzed using a very large and 
highly spectroscopically complete sample selected in the rest- 
frame NIR bolometric flux in the GOODS-N. They found that 
applying extinction corrections is critical when analyzing galaxy 
colors; nearly all of the galaxies in the green valley are 24pm 
sources, but after correcting for extinction, the bulk of the 24pm 
sources lie in the blue cloud. This correction introduces an av- 
erage shift in color of ~ 0.2 mag for the most extincted/star- 
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forming galaxies. However, to be consistent with the colors of 
the background galaxies no correction for intrinsic extinction is 
considered. 

AGN host-galaxies belong t o the red-sequenc e if their (U — 
V) color is above the threshold dBorch et al]l2006l) : 

(U - VWst-fmme > 0.277Log M, - 0.352 2 - 0.39 (5) 

Sources in the green-valley are approximately defined shifting 
this relation by 0.25 downward towards bluer colors. With these 
definitions, ~ 42% (108/255) and ~ 25% (63/255) of the to- 
tal sample are included in the red-cloud and the green-valley, 
respectively. For all sources the Specific Star-Formation Rate 
(SSFR) is estimated, defined as the ratio of the SFR per unit of 
galaxy stellar mass (SSFR=SFR/M«). The inverse of the SSFR, 
SS FR , is called "growth time" 



0.045SZ SO. 62 0.62< z SO. 88 0.88< z S3. 524 



SSFR 1 = MJM t 



(6) 



and corresponds to the time required for the galaxy to double its 
stellar mass, assuming its SFR remained constant. Actively star- 
forming galaxies are defined as sources with growth time smaller 
than the age of the Universe at their redshift (SSFR 1 < fHubbie), 
while sources with SSFR 1 larger than the age of the Universe 
can be considered passive galaxies (see also lFontana et al.f 2009: 
iBrusa et alJ2009h . Figure[T3"lshows SS FR' 1 as a function of the 
stellar mass in three different redshift bins for the AGN host- 
galaxies in the red-sequence, in the green-valley and in the blue- 
coud and for the zCOSMOS galaxies in same redshift ranges. 
The horizontal lines mark the age of the Universe at the two 
redshift boundaries of the chosen intervals. At face value, al- 
most all the sources in the red-sequence have SSFR' 1 larger 
than the age of the Universe at their redshift, which is consis- 
tent with passive galaxies. However, the value of SSFR' 1 has 
to be considered only as an approximate indication of the star- 
formation activity; in fact, there is some possible evidence of 
some residual star-formation, in red-cloud AGN host-galaxies, 
as witnessed by their morphologies. In the red-sequence 8 and 
28 sources are classified as ellipticals and SOs, respectively; all 
together they represent 34% of the host-galaxy population in 
the red-sequence. About 42% is represented by disk galaxies 
(both bulge-dominated and intemediate-bulge), which are prob- 
ably still forming stars but not at high rates. In fact, 15 over 108 
sources (~ 14%) have a detection at 70/vm and 5 have also a 
detection at 160yt/m (~ 6%). 

For these objects, the SFR inferred from the far-infrared de- 
tections is significantly higher than the SFR derived from the 
SED-fitting procedure. Indeed, an SED-fitting over the UV, op- 
tical and near-infrared bands is not always able to discriminate 
between the red continua of passive galaxies and those of dusty 
star-forming galaxies. Therefore, we decided to include another 
indicator in the present analysis bro adly following the procedure 
described in lCardamone et all (1201 Ol i.e., the (U -V) -(V - J) 
color diagram). Near-infrared emission can distinguish between 
red-passive or dust-obscured galaxies: given a similar 0.5/im 
flux, a star-forming galaxy has more emission near ~ 1/j.m than 
a passive galaxy. A sub-sample of galaxies is selected i n the 
same redshift range explored by ICardamone et al.l d2010l) . we 
find 92 AGN host-galaxies with 0.8 < z < 1.2. Fig. [141 shows 
both inactive galaxies and AGN host-galaxies in the same red- 
shift range and the thresholds considered to divide galaxies in 
the red-sequence and in the green-valley (we consider an aver- 
age redshift of 1 to define the threshold for the red-sequence 
and the green- valley). Thirty-five out of 92 AGN hosts are found 
to lie in the red-sequence (~ 38%) and 23 in the green-valley 
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Fig. 13. Inverse of the SSFR rate as a function of the stellar mass 
of the AGN host-galaxies in three different redshift bins for the 
zCOSMOS galaxies and for the Type-2 AGN sample in the red- 
sequence (red crosses), in the green-valley (green triangles) and 
in the blue-cloud (blue open circles). The horizontal lines mark 
the age of the Universe at the two redshift boundaries of the cho- 
sen intervals. 



(~ 25%); while for inactive galaxies about 32% and 21% lie in 
the red-sequence and green-valley, respectively. In Figure [T5lthe 
(U-V)-(V-J) color diagram for the 92 Type-2 AGN hosts is 
presented. From a preliminary analysis of the rest frame (U — V) 
against th e rest-frame (V — J) c olor (see Fig. Q3] but see also 
Fig. 2 in lCardamone et al.ll2010h , only ~ 9% of the AGN host- 
galaxies in the red-sequence and ~ 30% of AGN host-galaxies in 
the green-valley are moved in the region populated by dusty star- 
forming galaxies in the color-color diagram. To be compared 
with 20% AGN host-galaxies in the red-sequence and 75% AGN 
host-galaxies in the green-valley found by Cardamone and col- 
laborators. The fractions of dust-obscured galaxies among the 
red-cloud and green-valley A GN in our sample, a t 0.8 < z < 1 .2, 
are lower than those in the ICardamone et al.1 (1201 Ot) sample. 
However, the global fractions of AGN hosts, tentatively asso- 
ciated to passive galaxies, are very similar (~ 50%) in the two 
samples. The fractions for both AGN host-galaxies and inactive 
galaxies are reported in Table [2] 



9. Summary and Conclusions 

A detailed analysis of the SEDs of 255 spectroscopically identi- 
fied hard X-ray selected Type-2 AGN from the XMM-COSMOS 
survey is presented. In obscured AGN, the optical-UV nuclear 
luminosity is intercepted along the line of sight by the dusty 
torus and reprocessed in the infrared, so what we see in the 
optical-UV is mostly the light from the host-galaxy. On the one 
hand, this allows us to study the galaxy properties, on the other 
hand it makes difficult to estimate the nuclear bolometric power. 
An SED-fitting code has been developed with the main purpose 
of disentagling the various contributions (starburst, AGN, host- 
galaxy emission) in the observed SEDs using a standard^ 2 min- 
imization procedure (the starburst component is only used in the 
case of detection at 70//m). The code is based on a large set of 
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Table 2. AGN hosts and galaxies properties. 



Sample 


N 


Red- sequence 


Green-valley 


Blue-cloud 


0.045 < z < 3.452 


Type-2 AGN 


255 


108 (42%) 


96 (89%) P 
12(11%)D 


63 (25%) 


52 (82%) P 
11 (18%) D 


84 (33%) 


Galaxies 


8742 


2535 (29%) 


2306 (91%) P 
229 (9%) D 


1923 (22%) 


1596 (83%) P 
327 (17%) D 


4284 (49%) 


0.8 <z< 1.2 


Type-2 AGN 


92 


35 (38%) 


32 (91%) P 
3 (9%) D 


23 (25%) 


16 (70%) P 
6 (30%) D 


34 (37%) 


Galaxies 


1836 


587 (32%) 


569 (97%) P 
18 (3%) D 


385 (21%) 


269 (70%) P 
116 (30%) D 


864 (47%) 



Note - P=Passive, D=Dusty. 



0.8 Sz SIS 




Fig. 14. Distribution of the stellar masses as a function of the 
rest-frame (U - V) colors in the redshift range 0.8 < z < 1.2. 
The red dashed lin e represents the red sequence cut defined by 
iBorch et"aT] ((2006), while the green short dashed line defines an 
approximate green valley region, both lines are calculated at red- 
shift ~ 1. The points are color coded as in Fig. [12] 



starbu rst templates from ICharv & Elbazl ( 2001 ) and Dale et al] 
(l200lh . and galaxy templates from the lBruzual & Charlotl(l2003h 
code for s pectral synthesis m odels, while AGN templates are 
taken from lSilva et al] {2004). These templates represent a wide 
range of SED shapes and luminosities and are widely used in 
the literature. The total (nuclear) AGN bolometric luminosities 
are then estimated by adding the X-ray luminosities integrated 
over the 0.5-100 keV energy range to the infrared luminosity be- 
tween 1 and lOOOjum. The total X-ray luminosity is computed 
integrating the X-ray SED using the de-absorbed soft and hard 
X-ray luminosities. The SED is extrapolated to higher energies 
using the observed X-ray slope, and introducing an exponential 
cut-off at 200 keV. The total infrared luminosity is evaluated in- 




1 1.5 
rest-frame (V-J) 



Fig. 15. Distribution of Type-2 AGN hosts in the rest-frame 
(U - V) against the rest-frame (V - J) color. Color coded as 
in Fig. [13] Sources with the same best-fit galaxy template and 
the same extinction lie in the same position in the color-color 
diagram. Point size is keyed to the number of objects. 



tegrating the infrared AGN best-fit and then converted into the 
nuclear accretion disk luminosity applying the appropriate cor- 
rection factors to account for the geometry and the anisotropy of 
the torus emission. The reprocessed IR emission is considered 
to be a good proxy of the intrinsic disk emission and this is sup- 
ported by previous investigation s dPozzi et al.ll2007t Gandhi et 
al. 2009; Vas udevan et al. I l2oloh . In the distribution of the ratio 

r = Log (l 12 .3 ^obs/^12.3 ^.predicted \ see Eq. U the majority of 
the objects are within 2cr of the r distribution. The tail outside 
2cr and extending to high r includes 73 sources (with r > 0.5) 
for which the predicted mid-infrared luminosity is significantly 
lower than observed. We call "low-r" AGN all sources within 
2cr of the r distribution, while the "high-r" AGN sample is rep- 
resented by the sources deviating more than 2 <x. 
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Our main results are the following: 

1. The average observed SED is characterized by a flat X-ray 
slope, (F = 1.12), as expected for obscured AGN (not cor- 
rected for absorption), while in the optical-UV the observed 
light appears to be consistent with the host-galaxy emission. 
The average SED in the mid-infrared is more likely a combi- 
nation of dust emission from star-forming region and AGN 
emission reprocessed by the dust. 

2. The full sample is split into four bins of different X-ray and 
infrared luminosities and redshift. The shapes of the average 
SEDs in the optical bands are approximately the same in all 
luminosity and redshift bins. There is a stronger host-galaxy 
contribution at lower luminosity/redshift bins, where the av- 
erage SEDs have a typical galaxy shape. Moreover, there is a 
trend between X-ray and mid-infrared luminosity: the con- 
tribution of the AGN in the infrared (around 8 - 15//m) is 
higher at higher X-ray luminosities. 

3. Type-2 AGN appear to have smaller bolometric corrections 
than Type-1 AGN. At the same hard X-ray luminosity, 
43.30 < Log L[2-io]keV ^ 44.30, where both samples are well 
represented, we find that the median bolometric correction 
for Type-2 AGN (134 objects) is (k bol ) ~ 13 ± 1, to be com- 
pared with a median bolometric correction (k\, \) ~ 23 + 1 for 
Type-1 AGN (167 objects). The two averages are statistically 
different at the ~ 7 cr level. 

4. A clear separation in bolometric corrections for the low-r 
and the high-r samples is found. The relation provided by 
Gandhi and collaborators is valid for the majority of objects, 
while for 30% of the sample SED-fitting procedure may un- 
derestimate the non-nuclear contribution. At a given hard X- 
ray luminosity (43 < Log Lp-iojkeV ^ 44) the low-r sample 
has a median bolometric correction of (feboi) ~ H ± 1 (HO 
objects), to be compared with a median bolometric correc- 
tion for the high-r sample of (fcboi) ~ 26 + 3 (44 objects). 
The two median values for fcboi are statistically different at 
the ~ 5cr level. 

5. Host-galaxies morphologies and the stellar masses indicate 
that Type-2 AGN are preferentially hosted in galaxies which 
have a bulge, irrespective of the strength of the bulge or if the 
galaxy is on the red sequence or blue cloud, and with stellar 
masses greater than 10 10 M o . 

6. Almost all the sources in the red-sequence have SSFR 1 
larger than the age of the Universe at their redshift, which 
is consistent with passive galaxies. Following the same ap- 
proach as in Cardamone and collaborators (i.e., combining 
the rest-frame (U - V) vs Log M* and the rest-frame (U - V) 
vs (V - J) color diagrams), we find that, consistent with their 
results, ~ 50% of AGN hosts lie in the passive region of this 
diagram. In contrast from Cardamone et al. (2010), only ~ 
30% of AGN host-galaxies in the green-valley in our sample 
are consistent with dust-obscured sources in 0.8 < z < 1.2. 

It is clear that the mid and far-infrared parts of the SED are 
under-sampled with respect to the optical part. The ongoing 
Herschel survey over various fields at different depths (100//m 
and 160/im in the COSMOS field) and the upcoming ALMA sur- 
veys will allow us to gain an optimal multiwavelength coverage 
also in the far-infrared. 
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